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ABSTRACT 

The anisotropy of cosmic rays (CRs) in the solar vicinity is generally at- 
tributed to the CR streaming due to the discrete distribution of CR sources or 
local magnetic field modulation. Recently, the two dimensional large scale CR 
anisotropy has been measured by many experiments in TeV-PeV energy range 
in both hemispheres. The tail-in excess along the tangential direction of the 
local spiral arm and the loss cone deficit pointing to the north Galactic pole 
direction agree with what have been obtained in tens to hundreds of GeV. The 
persistence of the two large scale anisotropy structures in such a wide range of 
energy suggests that the anisotropy might be due to a global streaming of the 
Galactic CRs (GCRs). This work tries to extend the observed CR anisotropy 
picture from solar system to the whole galaxy. In such a case, we can find a new 
interesting signature, a loop of GCR streaming, of the GCR propagation. We 
further calculate the overall GCR streaming induced magnetic field, and find a 
qualitative consistence with the observed structure of the halo magnetic field. 

Subject headings: cosmic rays — diffusion — Galaxy: halo — ISM: magnetic 
fields 



1. INTRODUCTION 

Galactic magnetic field (GMF), Galactic cosmic ray (GCR) and the ordinary mat- 
ter are the basic components of the interstellar medium (ISM). These three constituents 
have comparable pressure and are bounded together by the electromagnetic force. The 
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GMF and the propagation of GCRs help to support the ordinary matter against the self- 
gravity. Conversely, the weight of the ordinary matter confines GMF and then GCRs in 
the Galaxy (Ferriere 2001). The dynamic balance process between these three constituents 
could give birth to new molecular-cloud complexes and ultimately trigger star formation 
(Mouschovias et al. 1974; Elmegreen 1982). The physical nature of and the interactions be- 
tween them are studied for decades. However, the fundamental question, the origin of them, 
is still open. Nowadays, due to the development of the observation technology and methods, 
the information about these constituents is more adequate, and we should approach a better 
understanding of this question. 

The GMF is commonly believed to be produced by a dynamo process with a seed field 
(e.g., see review (Brandenburg & Subramanian 2005) and references therein). Based on the 
measurements of the Faraday rotation of the linearly polarized radiation from pulsars and 
extragalactic radio sources (Simard-Normandin et al. 1981; Condon et al. 1998; Brown et al. 
2003; Han 2009; Taylor et al. 2009; Law et al. 2011), two kinds of structures are found in the 
halo magnetic field: the poloidal fields of dipole structure perpendicular to the plane in the 
Galactic center (GC) region and the toroidal fields of opposite directions above and below 
the Galactic plane (Han et al. 1997, 1999). The "AO" dynamo model (Han et al. 1997) could 
produce such a large scale magnetic field. Nevertheless, the origin of the GMF is still an 
open question. 

The GCR streaming, i.e. GCR propagation, could have contributions to the GMF. 
Parker (1992); Hanasz et al. (2009) suggested that the GMF could be generated by the dy- 
namo process driven by GCR streaming. Dolginov & Toptygin (2004) proposed a mechanism 
to produce the large-scale GMF by the electric current induced by GCRs. In their model, 
the current was theoretically calculated based on the diffusion law, and the result depends 
strongly on the assumption of the diffusion coefficient tensor. However, there is no direct 
measurement of the diffusion coefficient. It is better to use a direct method to estimate the 
GCR streaming. 

The anisotropy of GCRs should be a precise tool to probe the GCR streaming. Early in 
1935, Compton & Getting (1935) suggested that the relative movement between the GCR 
plasma and the observer due to the Galactic rotation could lead to a dipole anisotropy, known 
as the Compton- Getting (CG) effect. According to this effect, the direction and velocity of 
the CR streaming can be estimated by the phase and amplitude of the anisotropy respectively. 
Amenomori et al. (2006) did precise measurement above 300 TeV and concluded that the 
GCR plasma might corotate with the solar environment around the Galactic center with a 
velocity ^^220 kms~^. In such a manner, the CR streaming in the rest frame of the Galaxy 
could be deduced from the high-precision two-dimensional anisotropy map, which is observed 



- 3 - 



by many experiments in a wide energy range. 

In this work, we attempt to extend the anisotropy pattern observed in the solar vicinity 
to the whole Galaxy. The GCR streaming related to the anisotropy should exist in Galactic 
scale, forming a new picture of the GCR propagation. Furthermore, the global streaming 
of the GCR can be used to explore the contribution of GCR to the GMF. This Letter is 
organized as follows: In Section 2, we introduce the observational anisotropy and the GCR 
streaming in the extended picture. In Section 3, we estimate the contribution of the GCR 
to the GMF in the halo. In Section 4, we present discussions and the conclusions of this 
extension. 



2. ANISOTROPY AND COSMIC RAY PROPAGATION 

There is a long history to study the GCR anisotropy, with both underground /i de- 
tectors and ground-based air shower arrays (Aglietta et al. 1996; Munakata et al. 1997; 
Ambrosio et al. 2003; Antoni et al. 2004; Amenomori et al. 2005). Owing to the long-term 
observations and newly developed analysis methods, the high-precision two-dimensional map 
of the GCR anisotropy was obtained by Tibet air shower array (Amenomori et al. 2006), 
Super-Kamiokande (Guillian et al. 2007), MILAGRO (Abdo et al. 2009), ARGO (Zhang et al, 
2009) in the northern hemisphere, and IceCube (Abbasi et al. 2010) in the southern hemi- 
sphere. In this two-dimensional map, three broad structures are distinctly presented: an 
excess called the "tail- in", a deficit called "loss-cone" and an excess in the direction of the 
Cygnus region. 

Figure 1 (a) shows a schematic view of the GCR streamings in the solar system based 
on the northern hemisphere observations of the anisotropy. The three broad structures are 
shown as follows. 

The center of the "tail-in" component is at (Dec~-22°30', i?.A.~97°24'), close to the 
direction of the heliomagnetic tail {Dec^— 29°12' , R. A. ^90°2A') which is opposite to the 
proper motion direction of the solar system (Duldig 2001). This direction is quite close to 
the tangential direction of the local arm, indicating an inward GCR streaming along the 
local arm. A slight difference between the directions of the GCR streaming and the local 
arm might be due to the deflection by the heliomagnetic field. 

The center of the "loss-cone" deficit component points to the direction of the north 
Galactic pole (NGP). We interpret it as an outward streaming toward the NGP, the direction 
of which is perpendicular to the Galactic plane. In this picture, the GCR anisotropy of the 
southern hemisphere should show an excess toward the direction of the south Galactic pole 
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(SGP) (Z}ec~— 27°6', -R.y4.~12°48'). The observation from IceCube indeed showed such an 
excess (Abbasi et al. 2010). 

The excess component related to the Cygnus region peaks at (Dec~38°, 
i?.A.~309°). Similar with the above discussion, we expect this component to indicate an 
outward GCR streaming along the local spiral arm. 

The structures of the "tail-in" and "loss-cone" components are almost stable in the 
energy range between tens of GeV and hundreds of TeV, and insensitive to the solar activities 
(Amenomori et al. 2010). This indicates that the GCR anisotropy in this energy range may 
not be related to the heliospheric magnetic field, but of a large-scale origin. There is no 
reason that our solar location is special in the Galaxy. Therefore, we can extend the above- 
mentioned anisotropy pattern to the whole Galaxy. After that, the GCR plasma should have 
global streamings in three directions: inward and outward streamings along the spiral arms 
(The inward streaming related to the "tail-in" region is significantly larger than the outward 
streaming related to the Cygnus region. Thus, the overall streaming along the local arm is 
inward), and outward streaming perpendicular to the Galactic disk, as shown in Figure 1 



This global picture of the GCR propagation could be understood in the following way. 
The diffusion of GCRs is constrained by the magnetic field line, and the direction of the GCR 
streaming is decided by the gradient of the GCR intensity. Since most of the GCR sources 
are located in the inner Galactic disk, GCRs tend to diffuse out of the Galaxy in a direction 
perpendicular to the Galactic plane. There is a streaming back into the Galaxy from the 
direction of the spiral arm. Such a streaming may compensate the overall charge loss of the 
outward streaming, and makes the loop streaming of GCRs possible. Such a global loop 
streaming is not a new concept in astrophysics. For example, Alfven (1978) has discussed 
the electric double layers phenomenon and related loop electric current from laboratory scale 
to terrestrial, stellar and galactic scale. 

The streaming of the GCR plasma can form an electric current. The direction of the 
current is the same as the direction of the streaming because the GCR particles are mainly 
protons with positive charge (the fraction of relativistic electrons is less than 1%). The 
electric current density can be calculated by 



Here q is the mean charge of the GCR particles, we assume it to be the charge of proton in 
the following calculation for simplicity. E is the energy of the GCR particles. dN{E)/dE is 
the flux of the GCRs and a single power law spectrum (oc E~'^) with index 7 ^ 2.7 is used in 



(b). 
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this work. v{E) stands for the energy dependent streaming velocity of the GCRs, which can 
be estimated from the measurement of the anisotropy of GCRs. The integration is begun 
from 1 GeV and ended to 100 TeV as the contribution from higher energy is ignorable. 

The velocity of the GCR streaming is estimated by simply assuming the anisotropy of 
the GCRs is due to the CG effect (Compton & Getting 1935) 

^ « — (2) 

c 2 + 7 ^' ^ ' 

where c is the speed of light, and ^{E) is energy dependent anisotropy of GCRs (Parker 1992). 
According to the observations of various CR experiments, the amplitude of anisotropy in- 
creases with energy when the energy is below 1 TeV but nearly a constant when energy 
is between 1 and 100 TeV (Amenomori et al. 2005; Guillian et al. 2007). The anisotropy 
becomes to decrease when energy is higher than 100 TeV. For the purpose of an easy calcu- 
lation, the energy dependent amplitude is parameterized by the following expression: 

, , (axE'' (0.1 TeV <E <1 TeV) 
(a (1 TeV <E<100 TeV), 

where a ~ 0.0005, b ~ 0.44 are fitted from the observations. The energy is in a unit 
of TeV. For lack of the direct measurement, the amplitude of anisotropy below 0.1 TeV is 
estimated by the extrapolation of the above parameterization for energy between 0.1-1 TeV. 

As a result, the current density in the solar vicinity induced by the GCR plasma with 
energy above 1 GeV is ~ 3.1 x 10"^^ Am""^ . 

As shown in figure 1 (b), two kinds of Galactic electric currents are introduced: 

1) Inward currents along spiral arms 

The spiral arm of the Galaxy or the magnetic field in disk can be described by a loga- 
rithmic spiral function (Hou et al. 2009) 

R = re^^"'^. (4) 

Here {R, 9) are the polar coordinates and (p is the angle between the tangent and the radial 
line at the point (i?, 9). The complementary angle of if is called the pitch angle, and the mean 
pitch angle of Galactic spiral arm is 12° (Vallee 1995). When < < 37r and i? < 10 kpc 
the constant r is 1.35. 

We further assume that the inward current is mainly confined in the disk plane. It is ap- 
proximated by an exponential decay term in the Z direction (i.e., the direction perpendicular 
to the Galactic plane), 

j-(Z)=j7e-l^l/^\ (5) 
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where Hi is the characteristic decay length. Hi = lOOpc, the half-thickness of the inner disk 
(Ptuskin et al. 2009), is adopted in this work, is the inward current density in Galactic 
plane (Z=0). Using Equation (l)-(3), we already know the current density in solar vicinity, 
i.e. jcr{Z = lOpc), is 3.1 x 10"^^ Am""^, so is 3.4 x 10~^^ Am~'^ according to Equation 
(5). 

2) Outward currents perpendicular to the galactic disk 

Similarly, the outward currents to the NGP (for Z > 0) and the SGP (for Z < 0) are 
described by exponential functions in both the Z and R (radial) directions 

C\R, Z) = jOut^~\Z\/H,^^\R\/I^_ 

We adopt H2 = 4:kpc, and Rq = Akpc, in order to match the model of the halo magnetic 
field described by Sun & Reich (2010). It is reasonable for such a value of H2, due to 
the large uncertainty of the CR propagation halo size which could be larger than 10 kpc 
according to the present data (Putze et al. 2010). The exponential term of the radial direc- 
tion represents the feature that the density of CR sources is higher in the Galactic center 
region (Case & Bhattacharya 1996). j'q"* is the outward current density in GC, which is 
2.6 X 10~^^ Am~^ according to Equation (6), in the case that the current density in solar 
vicinity, j°"*(-R = 8.5kpc, Z = lOpc) is 3.1 x 10"^^ Am~^ (obtained by Equations (l)-(3), as 
mentioned above). 



3. GALACTIC HALO MAGNETIC FIELD 

To calculate the possible contribution of GCR to the GMF in the halo, firstly, we need 
to consider the electric currents from both the relativistic CRs and the thermal plasma 
(Dolginov & Toptygin 2004). The total current induced by these particles can be expressed 
as: 

J=3cr+Jl+Jth- (7) 

Here jcr is the current density induced by relativistic GCR particles, jf^ and j^i^ are the 
current densities induced by thermal protons and electrons respectively. 

As described in (Dolginov & Toptygin 2004), the diffusion coefficients for relativistic 
particles are not correlated with that for thermal particles. Thus, the two terms of the 
currents cannot cancel each other out, and should generate the magnetic field of the same 
order of magnitude. In this case, we can use the current generated by relativistic CR par- 
ticles as an order of magnitude approximation. It should be noted that the total current is 
divergence-free due to the charge conservation. 
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According to (Dolginov & Toptygin 2004), the evolution of the total GMF due to the 
electric currents of energetic particles can be solved by the following equation: 

^ = {i^m + T^turh)'^^B + V X (n X B) + i![^v X J,,. (8) 
ot c 

Here, Vm = (? jAixa is the local magnetic diffusion coefficient, z/.^ varies in different phases 
of the ISM: 3 x lO^^ cm^g-^ for cold (T - lOOiT), 3 x lO^o cm^s'^ for warm (T ~ lO^K) 
and 5 x 10^"^ cvr?s~^ for hot (T >~ 10^ Ji") phase (Dolginov & Toptygin 2004), and the ratio 
of filling factor for these phases roughly is 0.5 : 7.5 : 2 (Ferriere 2001). Therefore, the value 
averaged over the Galactic halo is ~ 3 x 10^*^ av?s~^ . u is the regular component of the 
plasma velocity. 

z/jtjrfe is the turbulent diffusion coefficient. It can be estimated using the formula Vturh ~ 
u'L/?)^ where L ^ lOOpc for typical component of the interstellar turbulence in the disk and 
u' ~ 10 kms~^ for the characteristic turbulent velocity. It gives fturb ~ 10^^ cim?s~^ in the 
disk. To the main concern of this work, i.e., the Galactic halo, vturb has been estimated to 
be 30 times higher by (Ruzmaikin et al. 1988), i.e., Vturb ~ 3 x 10^'^ cm^s~^. 

With these parameters, the diffusion time for the magnetic field to fill the Galaxy, i.e. 
the time to reach the steady state, can be estimated by r ~ / {uturb + J-'m) ~ 10^° years 
for Galactic radius R = lOkpc, which is in the same order as the age of the Galaxy. In 
other words, the current magnetic field of our galaxy is close to the stationary solution of 
the above equation, and the calculation of a steady magnetic field can be even simplified. 

In (Dolginov & Toptygin 2004), the above equation is obtained by applying an operation 
of V X on both side of the following equation: 

CE = -^^^jcr + {l^m + l^turb)^ X B -UX B. (9) 

c 

Considering that for the stationary solution, and the electric field associated with a plasma 
moving in a magnetic field is given by E = —u x B /c in the limit of large electrical con- 
ductivity, which is satisfied in most astrophysical circumstances (Parker 1979), Equation (9) 
can be rewritten as 

V X i? = - ^"^ 3.r. (10) 

C l^m ~r l^turb 

Comparing with the familiar Maxwell equation in the vacuum, we can see that the magnetic 
field in our case is equivalent to the solution for a largely suppressed electric current in the 
case of the vacuum. 



The result of the calculated magnetic field is shown in Figure 2. We can see clearly 
two components of the magnetic field, a poloidal component and a toroidal component. The 
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poloidal magnetic field is similar to a dipole pattern, with the direction from north to south, 
passing through the disk plane. The structure of the toroidal component is anti-symmetric, 
with a counter-clockwise direction in the north Galaxy and a clockwise direction in the south. 
Both structures are qualitatively consistent with the observations (Han et al. 1997, 1999). 

The radial profiles of the magnetic field at different Galactic heights (Z) are shown in 
Figure 3. We can see that the strength of the magnetic field decreases with the increase 
of Z. For the poloidal component, there is no observation yet to be compared directly. For 
the toroidal component, we give the parameterized description of observations, according to 
(Sun & Reich 2010). Our theoretical calculation shows a rough agreement of the structures 
with the observations. The magnetic field strength in our calculation is approximately 2% 
of the observation. 



4. DISCUSSION AND SUMMARY 

The observational anisotropy has a stable large-scale structure in the energy range 
from tens of GeV to hundreds of TeV. These structures imply three directions of the GCR 
streaming: inward and outward along the tangential direction of the spiral arm and the 
perpendicular direction to the Galactic pole. In this work, we extend the anisotropy observed 
locally to the whole Galaxy, which means the GCR streamings might exist at the global scale. 
A loop of the GCR streaming gives a new signature of the propagation of GCRs, and also 
provides contributions to the GMF. 

Qualitatively, the GCR streamings inferred from the anisotropy of the GCRs can gen- 
erate a large-scale magnetic field with poloidal and toroidal structures. These structures are 
consistent with the observations. This indicates that the extension of the local anisotropy to 
the whole Galaxy is probably reasonable. It will be helpful for further understanding of the 
origin of GCR anisotropy. Moreover, the anisotropy of the GCRs provides a new measure 
to study the Galactic electric current as well as a new window to understand the origin of 
the GMF. 

This analysis also indicates that the magnetic field in the Galactic halo may be partially 
contributed by the electric current induced by GCRs. The extension of the local anisotropy 
to the whole Galaxy is too simple in this work and the uncertainties introduced by the 
parameters used in the calculation are quite large, so the quantitative result only can be 
regarded as an order-of-magnitude estimate. 

We thank W. Y. Chen, S.F. Yang, and J. H. Xu for their questions, which stimulated 



- 9 - 



us to conduct this work. We are grateful to C. Q. Shen and Y. G. Xie for reading the 
manuscript and providing comments on the writing of the Letter. We express our gratitude 
to Q. Yuan, P. F. Yin, J. Zhang and M. M. Kang for the beneficial discussions. This work 
was supported by the Ministry of Science and Technology of China, the Natural Sciences 
Foundation of China (No. 10725524, 11105156) and by the Chinese Academy of Sciences 
(Nos. KJCX2-YW-N13, GJHZ1004). 

REFERENCES 

Abbasi, R. U., et al. 2010, ApJ, 718, 194 

Abdo, A. A., et al. 2009, ApJ, 698, 2121 

Aglietta, M., et al. 1996, ApJ, 470, 501 

Alfven, H. 1978, Astrophysics and Space Science, 54, 279 

Ambrosio, M., et al. 2003, Phys. Rev. D, 67, 042002 

Amenomori, M., et al. 2005, ApJ, 626, 29 

Amenomori, M., et al. 2006, Science, 314, 439 

Amenomori, M., et al. 2010, ApJ, 711, 119 

Antoni, T., et al. 2004, ApJ, 604, 687 

Brandenburg, A. & Subramanian, K. 2005, Physics Reports, 417, 1 
Brown, J. C, Taylor, A. R., & Jackel, B. J. 2003, ApJS, 145, 213 
Case, G. & Bhattacharya, D. 1996, A&AS, 120, 437 
Compton A.H. & Getting LA. 1935, Phys. Rev., 47, 817 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., Perley, R. A., Taylor, G. B., & 
Broderick, J. J. 1998, AJ, 115, 1693 

A. Z. Dolginov, & I. N. Toptygin 2004, A&A, 416, 411 

Duldig,M. L. 2001, PASA, 18, 12 

Elmegreen, B. G. 1982, ApJ, 253, 655 



- 10 - 

Ferriere, K. M. 2001, Reviews of Modern Physics, 73, 1031 
Guillian, G., et al. 2007, Phys. Rev. D, 75, 062003 

Han, J. L. 2009, in lAU Symp., 259, Cosmic Magnetic Fields:From Planets, to Stars and 
Galaxies (Cambridge: Cambridge Univ. Press), 455 

Han, J. L., Manchester, R. N., Berkhuijsen, E. M. & Beck, R. 1997, A&A, 322, 98 

Han, J. L., Manchester, R. N., & Qiao, G. J. 1999, MNRAS, 306, 371 

Hanasz, M., Woltanski, D., & Kowalik, K. 2009, ApJ, 706, 155 

Hon, L. G., Han, J. L., & Shi, W. B. 2009, A&A, 499, 473 

Law C. J. et al. 2011, ApJ, 728, 57 

Mouschovias, T. Ch., Shu, F. H., & Woodward, P. R. 1974, A&A, 33, 73 
Munakata, K., et al. 1997, Phys. Rev. D, 56, 23 

Parker, E. N., 1979, Cosmical Magnetic Fields: Their Origin and Their Activity, chap. 4(New 
York: Oxford Univ. Press) 

Parker, E. N. 1992, ApJ, 401, 137 

Ptuskin, V. S., Strelnikova, O. N. & Sveshnikova, L. G. 2009, Astrop. Phys., 31, 284 

Putze, A., Derome, L. & Maurin, D. 2010, A&A, 516, 66 

Ruzmaikin, A., Sokoloff, D. & Shukurov, A. 1988, Nature, 336, 341 

Simard-Normandin, M., Kronberg, P. P. & Button, S. 1981, ApJS, 45, 97 

Sun, X. H. & Reich, W. 2010, Research in Astronomy and Astrophysics, 10, 1287 

Taylor, A. R., Stil, J. M., & Sunstrum, C. 2009, ApJ, 702, 1230 

Vallee, J. P. 1995, ApJ, 454, 119 

Zhang J.L. et al. 2009, Proc. 31th Int. Cosmic Ray Conf. (Lodz) 
(http://icrc2009.uni.lodz.pl/proc/pdf/icrc0814.pdf) 



This preprint was prepared with the AAS IATJt;X macros v5.2. 



- 11 - 




Fig. 1. — (a) Schematic view of the anisotropy of the GCRs observed at the solar position. 
The three red arrows correspond to three anisotropy components, namely, the "tail-in" ex- 
cess, the "loss-cone" deficit, and the Cygnus excess, respectively. The directions of their 
centers are noted as Itaii^in, hoss-cone and I cygnus respectively, haii-in and I cygnus are approx- 
imately opposite along the local arm. Iioss-cone is perpendicular to the Galactic plane, (b) 
Cartoon picture of the OCR streamings after the extension of the local anisotropy to the 
whole Galaxy. Inward streamings (j*") are along the spiral arms and outward streamings 
(j""*) are perpendicular to the Galactic disk. 
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Fig. 2. — Structure of the Galactic halo magnetic field contributed by the GCRs. The 
GMF in the halo is composed of a poloidal magnetic field and a toroidal magnetic field. 
The direction of the poloidal field is from Galactic north to south, passing through the disk 
plane. The toroidal field has reverse directions below and above the Galactic plane. The 
arrows indicate the directions of the magnetic fields. 
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Fig. 3. — Distribution of the poloidal (a) and toroidal (b) magnetic field strength. The 
solid line is the strength of the magnetic field calculated by our model as a function of 
the Galactocentric radius (R) at a given height (Z). The dashed lines in (b) present the 
strength distribution of the toroidal magnetic field derived from the parameterized model 
(Sun & Reich 2010). The strength of the magnetic field (solid line) is the magnetic field 
calculated by our model multiplied by 50 for the convenience to compare the distribution 
behavior. The solid lines and the dashed lines from top to bottom correspond to the height 
at 2, 4, 6, and 8 kpc, respectively. 



